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ABSTRACT 

Using six high resolution dissipationlcss simulations with a varying box size in a flat 
LCDM universe, we study the mass and redshift dependence of dark matter halo 
shapes for Mvir — 9.0 x 10^^ — 2.0 x 10^^, over the redshift range z = — 3, and for 
two values of as — 0.75 and 0.9. Remarkably, we find that the redshift, mass, and 
(78 dependence of the mean smallest-to-largest axis ratio of halos is well described by 
the simple power-law relation (s) = (0.54 ± 0.02)(Afvir/M*)~°'°''°^°'°°^: where s is 
measured at 0.3i?vir and the z and as dependences are governed by the characteristic 
nonlinear mass, Af* — M^,{z,as). We find that the scatter about the mean s is well 
described by a Gaussian with a ^ 0.1, for all masses and redshifts. We compare our 
results to a variety of previous works on halo shapes and find that reported differences 
between studies are primarily explained by differences in their methodologies. We 
address the evolutionary aspects of individual halo shapes by following the shapes of 
the halos thr ough ~ 100 sn apshots in time. We determine the formation scalefactor Qc 
as defined bv lWechsler et a l. (200^ and find that it can be related to the halo shape 
at z = and its evolution over time. 

Key words: cosmology: theory — galaxies: formation — galaxies: halos — large-scale 
structure of universe 



1 INTRODUCTION 



A generic prediction of cold dark matter (CDM) theory is 
the process of bottom up halo formation, where large halos 
form from the mergers of smaller halos, which are in turn 
formed from even smaller halos. This is a violent process 
and it violates most of the assumptions that go into the 
spherical top-hat collapse model of halo formation which is 
often used to describe halos. Since mass accretion onto halos 
is often directional and tends to be clumpy, one would not 
expect halos to be spherical if the relaxation time of the 
halos were longer than the time between mergers and/or if 
the in-falling halos came along a preferential direction (such 
as along a filament). In both theoretical modelling of CDM 
and observations, halos are found to be very non-spherical. 
In fact, spherical halos are rare. Therefore, the analysis of 



halo shapes can give us another clue to the nature of the 
dark matter and the process of halo and galaxy formation. 

One way of quantifying the shape of a halo is to go 
one step beyond the spherical approximation and approx- 
imate halos by ellipsoids. Ellipsoids are characterised by 
three axes, a, b, c, with a > b > c, which are normally de- 
scribed in terms of ratios, s = c/a, q = b/a, and p = c/b. El- 
lipsoids can also be described in terms of three classes, which 
have corresponding ratio ranges: prolate (sausage shaped) 
ellipsoids have a > b ~ c leading to axial ratios of s ~ g < p, 
oblate (pancake shaped) ellipsoids have a ~ fe > c leading 
to axial ratios of s « p < 5, and triaxial ellipsoids are in be- 
tween prolate and oblate with a > b > c. Additionally, when 
talking about purely oblate ellipsoids, a = 6, it is common 
to use just q, since s and q are degenerate. 

There have been many theoretical papers pub- 
lished over the years which examined the subject 
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of halo 
eludes 



shapes. The earl y work on the subject in- 



Barnes fc Efstathioul dlOSTT): iDubinski fc Carlberd 
Katz ('1991*); 'Warren ct al. ('l992^;'Dubinski' ('l994|); 



IJing et al. (1995): Tormen (1997); Thoma s ot al. (1998). AU 
of these works agreed that h alos are ellipsoidal, but oth - 
erwise differ in several details. iDubinski fc Carlberd (|l99ll) 
found that halos have axial ratios of s ~ 0.5 in the interior 
and b eeome more spherical at larger radii, while iFVenk et all 
il988l) found that halos are slightly more spherical in the 
centres. iThomas et al.l (Il99^ claimed that larger mass ha- 
los have a slight tendency to be more spherical, where more 
recent simulation results find the opposite. Despite these dis- 
agreements many of the early authors g ive us clues into the 
nature of halos shapes. IWarren et al.l (Il992h and iTormenI 
1^93) showed that the angular momentum axis of a halo 
is well correlated with the smallest axis, c, although, as 
most of these authors pointed out, halos are not rotation- 
ally supported. This therefore has led many to conclude that 
the shapes are supported by anisotropic velocity dispersion. 
iTormenI il997l) took it one step further and found that the 
velocity anisotropy was in turn correlated with the infall 
anisotropy of merging satellites. Most authors found that 
the axial ratios of halos are ~ 0.5 ± 0.1 and that halos tend 
to be prolate as opposed to oblate in shape. The most likely 
source of the disagreement in these works and in the more 
recent works we describe below is the different methods used 
often coupled with inadequate resolution. 

More rec ent stud ies of hal os' shape w er e per - 
formed bv iBuUockl (I2p02h: Ijing fc Sutol J2002 " 

" " 3j); iBj 



ISgrinee^^^i^^^^^Muis^ (120041): iBailin fc Steinmet j 
^2005^ : lKasun fc EvrardI J2004) : iHopkins et alJ ^2005^ . The 

results of these authors differ, in some cases, considerably. 
One of the goals of this paper is to carefully examine the 
differences in the findings presented by the above authors. 

All of t he a forementioned publications (except 
ISpringel et al.l |2004) analyse simulations with either no 
baryons or with adiabatic hydrodynamics. This is both due 
to the cost associated with performing self consistent hy- 
drodynamical simulations of large volumes with high mass 
resolution and due to the fact that very few cosmological 
simulations yet produce realistic galaxies. Nonetheless, we 
know that the presence of baryons should have an effect 
on the shapes of halos due to their coUisional behaviour. 
Three recent papers have attemp t ed to examine the effect 
of barvons iSpringel et al.l 120041: iKazantzidis et al.l l2004t 
iBailin et al1l2005ft . In lSoringel et al.l (l2004^ the same simu- 
lations were done using no baryons, adiabatic baryons and 
baryons with cooling and star formation. In the first two 
cases there was very little difference, but with the presence 
of cooling and star formation the halos became more 
spherical. The radial dependence of shape also changes such 
that the halo is more spherical in the centre. At 7i > 0.3-Rvir 
the axial ratio s increases by < 0.09, but in the interior the 
increase As is as much as 0.2. In an independent study, 
IKazantzidis et al] (|2004|) found an even larger effect due 
to baryonic cooling in a set of 11 high resolution clusters. 
At R = 0.3J?vir the authors found that s can increase by 
0.2 — 0.3 in the presence of gas cooling. The extent of the 
over-cooling problem plaguing these simulations is still 
uncertain. This amount of change in the shape should be 
viewed as an upper l i mit. The most recent work on the 
subject is lBailin et al.l (l2005h . who concentrate more on the 



relative orientation of the galaxy formed at the centre of 
eight high resolution halos than on the relati ve sphericity 
of the halos. Despite this, from Figure 1 of iBailin et alJ 
^005) it seems that they would also predict an increase 
of ~ 0.2 for s. It is still useful to study shapes of halos 
without baryonic cooling. Cooling and star formation in 
simulations is still a very open question, making the effect 
of the cooling uncertain. We show in Paper II that our 
simulations without cooling match shapes of X-ray clusters. 

Measurements of the shapes of both cluster and 
galaxy mass halos through varied observational tech- 
niques are increasingly becoming available. There 
have been many studies of the X-ray morpholo- 
gies of clusters (see iMcMillan. Ko walski fc UlmeI^ 

[19891: iMohr. Evrard. Fabricant fc GeUeil Il99a 

'Kol okotronis. Basilakos. Plionis fc Georgantopouloj 1200 j) 
which can be directl y related to the shape of the inner par t 
of the cluster halo jLee fc Sutd 120031: iBuote fc Xulll997D . 
For a review of X-ray cluster shapes and the latest results, 
see Florcs ct al] j^QO^), hereafter referred to as Paper II. 

There has also been important new information on the 
shape of galaxy mass halos, in particular our own Milky Way 
halo. Oiling fc Merrificld (2000) concluded that the host 
halo around the Galactic disk is oblate with a short-to-long 
axial ratio of 0.7 < q < 0.9. Investigations of Sagittarius' 
tidal streams have led to the conclusion that the Milky Way 
halo is oblate and nearly spherical with q > 0.8 ('ibata et ajj 
^001; Maicwski ot al. 2003; Martmcz-Dclgado ct al. 2004^^^ 
However, by insp ecting M giants within the leading stream 
iHelmi (|200J) and iLaw. Johnston fc M^ewskJ i2005|) found 
a best fit prolate halo with s = 0.6. iMerrifieldl (|2004^ sum- 
marises the currently reliable observations for galaxy host 
halo shapes using multiple techniques and find that the ob- 
servations vary a lot. 

Another method for studying shapes of halos at 
higher redshift is galaxy-galaxy weak lensing studies. 
Analy sing data taken with the Canada- France-Hawaii tele- 
scope, iHoekstra. Yee fc Gladder^ i2004h find a signal at a 
99.5% confidence level for halo asphericity. They detect an 
average projected ellipticity of (e) = (1 — 920) ~ 0.201^q q5, 
corresponding to s = 0.66lQ og, for halos with an average 
mass of 8 X 10^^ Mq. Ongoing studies of galaxy-galaxy 
weak lensing promise rapidly improving statistics from large 
scale surveys like the Canada-France Legacy survey. 

This paper is organised as follows: In Section |21 we de- 
scribe the simulations, halo finding method, and halo prop- 
erty determination methods used in this study. In Section 
13 we discuss the method used to determine the shapes of 
halos. In Section 2] we examine the mean axial ratios from 
our simulations and their dependence on mass, redshift and 
(Tg. We then examine the dispersion of the axial ratio ver- 
sus mass relation. We briefiy discuss the shape of halos as 
a function of radius and then examine the relationship of 
the angular momentum and velocity anisotropy to the halo 
shape. In Sectionj^we examine the relationship between the 
formation history of halos and their present day shapes. In 
Section |S| we compare our results to those of previous au- 
thors and explain the sources of the differences. In section 
[7| we examine the observational tests and implications of 
our findings. Finally, Section |H| is devoted to summary and 
conclusions. 
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Table 1. Simulation parameters 



Name 






^box 




r?ip 


^pcak 




Af,(10i2 


h-l Mq) 














/i-l Mq 


kpc 


z = 


2 = 1 


2 = 2 


2 = 3 


L8O0.75 


0.75 


0.030 


80 


512^ 


3.16 X 10* 


1.2 


3.0 


0.11 


0.0046 


0.00027 


L80o.9a 


0.9 


0.045 


80 


512^ 


3.16 X 10* 


1.2 


8.0 


0.35 


0.019 


0.0013 


LSOo.gt, 


0.9 


0.045 


80 


512^ 


3.16 X 10* 


1.2 


8.0 


0.35 


0.019 


0.0013 


L2OO0.9 


0.9 


0.030 


200 


256^ 


3.98 X 10^0 


5.0 


8.6 


0.41 


0.023 


0.0018 


LI2O0.9 


0.9 


0.045 


120 


5123 


1.07 X 10^ 


1.8 


8.0 


0.35 


0.019 


0.0013 


L120o.9r 


0.9 


0.045 


40 sphere 


~ 256^ 


1.33 X 10* 


0.9 


8.0 


0.35 


0.019 


0.0013 



2 SIMULATIONS 

2.1 The Numerical Simulations 

All our simulations (see Table were performed with 
the Adaptive Refinement Tr e e (AR T) N-body code of 
iKravtsov. Klvpin fc Khokhlovl lll997l) which implements 
successive refinements in both the spacial grid and time step 
in high density environments. We analyse the shapes of ha- 
los and their merger histories in the concordance fiat ACDM 
cosmological model: Q.m = 0.3 = 1 — Q,h, h — 0.7, where 
Qm and Qa are the present-day matter and vacuum energy 
densities in units of critical density and h is the Hubble pa- 
rameter in units of 100km s~'^Mpc~^. The power spectra 
used to generate the initial conditions for the simulations 
were determined from a direct Boltzmann code calculation 
(courtesy of Wayne Hu) . 

To study the effects of the power spectrum normal- 
isation and resolution, we consider five simulations of 
the ACDM cosmology. The first simulation (L8O0.75) fol- 
lowed the evolution of 512^ = 1.34 x 10* particles in a 
80h~^ Mpc = 114.29Mpc box and was normalised to as — 
0.75, where cts is the rms fluctuation in spheres of 8h~^ Mpc 
comoving radius. The second simulation (L8Q0.9) is an ex- 
act replica of the L8O0.75 simulation with the same random 
number seed, but the power spectrum was normalised to 
as = 0.9. The flrst simulation was also used to study the halo 
oc cupation distribution a nd t he physics of galaxy cl ustering 
bv lKravtsov et alJ J2004 and lZentner et al.l ^M5^. Unfor- 
tunately, both of these simulations were generated with a 
power spectrum which had a little more than average power 
on large scales. This may happen when generating power 
spectra due to cosmic variance. The simulation is still a good 
representation of a volume in the Universe, but to avoid be- 
coming non- linear on large scales, the second simulation was 
stopped at z = 0.1. Due to the lower normalisation of the 
L8O0.75 box it was allowed to run until z = 0. We use these 
two simulations to study the effects of the spectrum normal- 
isation, but to achieve better statistics and make predictions 
for as = 0.9 at 2 = we also include another simulation of 
the same size and resolution (L80o.g6). The fourth simula- 
tion (L2OO0.9) followed the evolution of 256^ = 1.68 x lO'^ 
particles in a 200/i"^ Mpc = 285.7Mpc box. The flfth sim- 
ulation LI2O0.9 followed the evolution of 512^ particles in 
a 120/1"^ Mpc = 171.43Mpc box and was normalised to 



as = 0.9. This simulation is used for several purposes: firstly 
to achieve better statistics for rare high mass objects, and 
secondly as the basis for the sixth simulation. The sixth sim- 
ulation is a resimulated Lagrangian subregion of the LI2O0.9 
box corresponding at z = to a sphere in position space 
of diameter D — 40/i~^ Mpc. The initial conditions of the 
LI2O0.9 box initially contained 1024^ particles which were 
combined into 512^ particles used for the initial simulation. 
The Lagrangian subregion was then chosen and the original 
higher resolution particles of mass nip = 1.33 X 10*/i-i Mq 
within this region, corresponding to 1024^ particles in the 
box, were followed from the initial time step, Zi — 40. The 
high mass resolution region was surrounded by layers of par- 
ticles of increasing mass with a total of five particle species 
in order to preserve the large scale gravitational field. Only 
the regions containing the highest resolution particles were 
adaptively refined. The maximum level of refinement in the 
simulation corresponded to a peak formal spatial resolution 
of 0.9/i~^kpc. For more details about the multi-mass tech- 
nique consult 'Klvp in et alJ i200 j) . The subregion was cho- 
sen not to contain any halos above Mvir > W^^h~^ Mq in 
order to increase the statistics of isolated galaxy mass halos. 



2.2 Halo Identification and Classification 

A variant of the Bound Density Maximum (EDM) algo- 
rithm is used to identify h alos and subhalos in our sim- 
ulations telvpin et aLllTggsl) . The details of the algorithm 
a nd parameters being u sed in the halo finder can be found 
in lKravtsov et all (120041) . We briefiy describe the main steps 
in the halo finder here. First, all particles are assigned a den- 
sity using the smooth algorithm^, which uses a symmetric 
SPIf (Smoothed Particle Hydrodynamics) smoothing ker- 
nel on the 32 nearest neighbours. Density maxima are then 
identified which are separated by a minimum distance of 
rmin = 50h~^ kpc, defining the minimum distinguishable 
separation of halos and subhalos. Using the maxima as cen- 
tres, profiles in circular velocity and density are calculated in 
spherical bi ns. Unbound particle s are removed iteratively as 
described in lKlvpin et alJ 1^9^. The halo catalogue used is 

^ To calculate the density we use the publicly available 
code smooth; ht tp: //wgw-hpcc . astro .Washington. edu/tools7| 
tools.html 
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complete for halos with > 50 particles. This corresponds to 
a mass below which the cumulative mass and v elocity func- 
tions begin to flatten fsee lKravtsov et alJ lj20o3) for details) . 

The halo density profiles are constructed using only 
bound particles an d they are fit by an NFW profile 
jNavarro et al.lll996l) : 



where 



PNFw{r) 



{r/r,){l + r/rsY 



(1) 



where is the radius at which the log density profile has 
a slope of —2 and the density is Ps/4. One of the pa- 
rameters, or pg, can be replaced by a virial parameter 
(-RvirjMvir, Or Vvii) defined such that the mean density in- 
side the virial radius is Avir times the mean universal density 
po(z) — flm{z)pc{z) at that redshift: 



47r 3 



where pc{z) is the critical density, and 



I8n^ +82{n„ 



39{n„,{z) 



^m{z) 



(2) 



(3) 



from iBrvan fc NormanI Jl998l) with Avir(O) « 337 for the 
ACDM cosmology assumed here. The NFW density profile 
fitting is performed using a minimisation algorithm. The 
profiles are binned logarithmically from twice the resolution 
length (see Tabled out to iisoo, the radius within which 
the average density is equal to 500 times the critical density 
of the universe. The choice of this outer radius is motivated 
by [Tasitsiomi et al. (2004) who showed that halos are well 
relaxed within this radius. The binning begins with 10 ra- 
dial bins. The number of bins is then reduced if any bin 
contains fewer than 10 particles or is radially smaller than 
the resolution length. This reduction of bins is continued un- 
til both criteria are met. Fits using this method have been 
compared to fits determined using different merit functions, 
such as the maximum deviation from the fit as described in 
iTasitsiomi et all ll2004ll and it was found that they give very 
similar results for individual halos. After fitting the halos the 
host halo and subhalo relationship is determined very sim- 
ply. If a halo's centre is contained within the virial radius 
of a more massive halo, that halo is considered a subhalo 
of the larger halo. All halo properties reported here are for 
halos which are determined to be isolated or host halos (i.e., 
not subhalos). 



3 METHODS OF DETERMINING SHAPES 

There are many different methods to determine shapes of 
halos. All methods model halos as ellipsoidal with the eigen- 
vectors of some form of the inertia tensor corresponding to 
the axes c<b<a{s = c/a and q = b/a). The two forms of 
the inertia tensor used in the literature to determine shape 
are the unweighted. 



n 

and the weighted (or reduced). 



(4) 



(5) 



\Jx'i + yl/q'^ + zl/s'^, 



(6) 



is the elliptical distance in the eigenvector coordinate system 
from the centre to the nth particle. In both cases the eigen- 
values {\a < Ai, < Ac) determine the axial ratios described 
at the beginning of SectionQwith (a, b, c) = y/Xa, A;,, Ac. The 
orientation of the halo is determined by the corresponding 
eigenvectors. 

One would like to recover the shape of an isodensity sur- 
face. The method used here begins by determining / with 
s = 1 and q = 1, including all particles within some radius. 
Subsequently, new values for s and q are determined and 
the volume of analysis is deformed along the eigenvectors 
in proportion to the eigenvalues. There are two options to 
choose from when deforming the volume. The volume within 
the ellipsoid can be kept constant, or one of the eigenvectors 
can be kept equal to the original radius of the spherical vol- 
ume. In our analysis of shapes, the longest axis is kept equal 
to the original spherical radius. After the deformation of the 
original spherical region, / is calculated once again, but now 
using the newly determined s and q and only including the 
particles found in the new ellipsoidal region. The iterative 
process is repeated until convergence is achieved. Conver- 
gence is achieved when the variance in both axial ratios, s 
and q, is less than a given tolerance. 

The analysis presented here begins with a sphere oi R = 
0.3i?vir, and keeps the largest axis fixed at this radius unless 
otherwise stated. For determining halo shapes accurately we 
limit our analysis to isolated halos with Np > 7000 within 
7?vir. This corresponds to Mvir > 2.21 x lO^^/i"^ M© for 
the 80ft" ^Mpc box simulations, Mvir > 7.49 x 10^^ Mq 
for the 120/i~'^Mpc box simulation, and Mvir > 9.3 x 
W^^h~^ Mq in the resimulated region of the 120ft"^Mpc 
box. For a discussion of our resolution tests, see Appendix 

El 



4 SHAPES AS A FUNCTION OF HALO MASS 

The simulations analysed here enable us to analyse halos 
spanning a mass range from galaxy to cluster sized objects. 
These data provide an opportunity to study the variation 
of shape and its intrinsic scatter with halo mass. Various 
statistics are used to derive robust estimates of the depen- 
dence of shape on halo-centric radius. Combining the de- 
tailed spatial and dynamical information from the simula- 
tions we can relate quantities like angular momentum or 
velocity anisotropy tensor to the shape and the orientation 
of the halo. In this section we aim to present a comprehen- 
sive analysis of the properties of halos at all redshifts. In 
section 13 we will address evolutionary aspects of individual 
halo shapes. 



4.1 Median Relationships for Distinct Halos 

We begin by fitting the mass dependence of halo shape and 
find that the mean value of the axial ratio s = c/a decreases 
monotonically with increasing halo mass as illustrated in 
Figure Q In other words, less massive halos have a more 
spherical mean shape than more massive halos. Since we use 
four different simulations (L8O0.96, LI2O0.9, L120o.9r, and 
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■ 10'2 1013 iQii iQis 

(Mo h-') 

Figure 1. Mean axial ratios s = c/a for four simulations of 
different mass resolution are presented with a fit (solid black 
lines) given by Equation J7J and dispersion of 0.1. The trian- 
gles, squares, solid circles and X symbols are the average s for 
a given mass bin. The solid circles have been shifted by 0.05 in 
log for clarity. The open circles and error bars are the best fit 
Kolmogorov-Smirnov mean and 68% confidence level assuming a 
Gaussian parent distribution. The dashed lines connect the raw 
dispersion for each point and the coloured solid lines are the best 
fit (KS test) dispersion. (See the electronic edition for colour ver- 
sion of the figure) 

L2OO0.9) with varying mass and length scales we are able to 
determine (s)(Mvir) over a wide mass range. We find that 
over the accessible mass range the variation of shape with 
halo mass is well described by 

{s)(Mvi.,^ = 0) = a(^^y (7) 

with best fit values 

Q = 0.54 ± 0.03, P = -0.050 ± 0.003. (8) 

The parameters, a and /3 were determined by weighted 
minimisation on the best fit mean data points determined 
via Kolmogorov-Smirnov (KS) analysis assuming a Gaussian 
distribution within a given mass bin (see section 1^^ . M*(z) 
is the characteristic nonlinear mass at z such that the rms 
top-hat smoothed overdensity at scale a{M,,z) is 5c = 1.68. 
The M, for 2 = is 8.0 x 10^^/i~^ Mq for the simulations 
with ~ 0.045 and 8.6 x 10^^ for the simulations with 
flfc — 0.03. Only bins containing halos above our previously 
stated lower bound resolution limit were used and only mass 
bins with at least 20 halos were included in the fit. This 
work extends the mass ra nge of the similar relationships 
found by previo us authors il.Ting fc Sut J2OO2I: lBullockll2002l: 
ISnringel et al.l 2004; Kasun fc Evrardl2004! ) : we compare our 
results with these previous works in Section |S| 




K^r (Mo h-') 

Figure 2. {s){M) for z = 0.0,1.0,2.0,3.0. The binning is the 
same as in Figure but now for many different redshifts. The 
solid line is the power-law relation set out in Equation Q. 
The LI2O0.9 points are shifted by 0.05 in log for clarity. The 
Springel data agrees quite well with our data and model for 
z = 0.0, 1.0,2.0. 

4.2 Shapes of Halos at Higher Redshifts 

The use of M« in the Equation O alludes to the evolution 
of the (s)(Mvir) relation. After examining the (s)(Mvir) re- 
lation at higher redshifts, we find that the relation between 
(s) and Mvir is successfully described by Equation O with 
the appropriate Mt{z). The M* for z = 1.0, 2.0, and 3.0 
are 3.5 x lO", 1.8 x 10^", and 1.3 x lO'^h'^ Mq respectively 
for the simulations with Qi, = 0.045. We present our results 
for various redshifts in Figure |5| from the L120o.9r, L8O0.9, 
LI2O0.9 and Z/200o.9 simulations. We have also included data 
points provided to us by Springel (private communication) 
in Figure |5| for comparison, which from a more complete 
sample than the data presented in Spri ngel ct al. . (.2004!) and 
are for shapes measured at 0.4J?vir. 



4.3 Dependence on erg 

Of the parameters in the ACDM cosmological model the pa- 
rameter which is the least constrained and the most uncer- 
tain is the normalisation of the fluctuation spectrum, usually 
specified by as- Therefore, it is of interest to understand 
the dependence of the (s)(Mvir) relation on as- Since Af, 
is dependent on as the scaling with M« in Equation Q 
may already be sufficient to account for the as dependence. 
As stated in Section |51 L8O0.75 and L80o.9a were produced 
with the same Gaussian random field but different values 
for the normalisation. Therefore the differences between the 
two simulations can only be a result of the different values 
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Figure 3. (s) vs Al with different values of erg. Different values of 
as predict different values for the (s) vs M relationship. Here one 
can see that a universe with a lower erg produces halos which are 
more elongated, although the power-law relationship (Equation 
171 1 remains valid, as shown by the agreement between the points 
and the lines representing this prediction. 



for as. As Figurel^iUustrates, the two simulations do indeed 
produce different relations. We find that the M* dependence 
in Equation ||7J is sufficient to describe the differences be- 
tween simulations of different as- One should expect this 
from the result of the previous subsection, that the redshift 
evolution was also well described by the M, dependence. 
The values of M. for z = 0.1 are 5.99 x 10^^ for as = 0.9 
and 2.22 x 10^^ for erg = 0.75. The value of M* for as = 0.75 
at 2 = 1 and 2 are 1.09 x 10^^ and 4.57 x 10^ respectively. 
A simple fit to the redshift dependence of A/* in these cos- 
mologies is log(M.) = A - Blog(l + z) - C(log(l + z))^, 
with A{B, C) = 12.9(2.68, 5.96) for erg = 0.9 and A{B, C) = 
12.5(2.94, 6.28) for erg = 0.75, and accurate to within 1.6% 
and 3.1%, respectively, for z < 3. 



4.4 Mean - Dispersion Relationship 

In the previous subsections we used the best KS test fit 
mean, assuming a Gaussian parent distribution, as an esti- 
mate of the true mean of axial ratios within a given mass 
bin. In this subsection we examine the validity of this as- 
sumption, and test wh ether the disp e rsion has the mass de- 
pendence suggested by Ijing fc Sutd ||2002|) . In Figure 2] we 
present the distribution of s in the six bins from FigureQfor 
LI2O0.9. In each of the plots we have also included the KS 
best fit Gaussian, from which the mean was used to deter- 
mine the best fit power-law in Equation Q. The error-bars 
on the mean indicated in Figure are the 68% confidence 
limits of the KS probability. The limits are determined by 



Figure 4. The distribution of s in the L120o.9i, simulation in the 
mass ranges indicated. The number of halos in each bin is also 
indicated. The Gaussian fit shown for each graph is the best fit, 
based on a KS test analysis. 



varying the mean of the parent distribution until the KS 
probability drops below 16% for greater and less than the 
best fit value for the dispersion in each mass bin. The values 
from this analysis corresponding to the distributions in Fig- 
ure 2] can be found in Table |21 In Figure 0] the lowest mass 
bin, which also contains the most halos, is well fit by a Gaus- 
sian. This is seen in Table|5|not only by the best fit KS prob- 
ability, but also by the small range of the confidence limits. 
The higher mass bins are consistent with having Gaussian 
parent distributions though the parent distributions' values 
for the mean and dispersion are not as well constrained. 
There is no indication of a structured tail to lower values of 
s, but Table 121 indicates that the distributions have negative 
skewness. This arises from a small number of halos with very 
low values of s, which are always determined to be ongoing 
major mergers with very close cores. 

ping & Suto (2002) found that the distribution of s 
within a given mass bin is Gaussian. They found no indi- 
cation of a tail or any low values of s. This is most likely 
due to t heir treatment of halos with multiple cores (see sec- 
tion ISll. lBullockl (|20o3) found a large tail to low values of s 
using R = Rvii- After repeating our analysis at _R = J?vir we 
find the exact opposite. We find even less indication of a tail 
than in the distributions shown in Figure |1] The difference 
is most likely due to the ce ntres of halos determ i ned b y the 
different halo finders used. iBailin fc Steinmet j (j200^ find 
a more subtle but significant tail to lower values of s. This 
is most likely just a side effect of combining all mass bins 
into one histogram. If the histogram were divided into bins 
over smaller ranges in mass, this tail would be seen as a con- 
sequence of the combination of Gaussian distributions with 
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Table 2. Kolmogorov-Smirnov Best Fit Values 



Mass 
(h-^ Mq) 










era 




KS Prob 


4.8 X 10^2 


0.583 (H 


-0.003 


-0.003) 


0.108 (H 


-0.006 


-0.005) 


0.80 


9.6 X 1012 


0.554 (H 


-0.006 


-0.006) 


0.110 (H 


-0.005 


-0.007) 


0.61 


1.92 X 10^^ 


0.518 (H 


-0.009 


-0.004) 


0.094 (H 


-0.007 


-0.013) 


0.72 


3.84 X 10^^ 


0.519 (H 


-0.014 


-0.013) 


0.108 (H 


-0.016 


-0.030) 


0.95 


7.68 X 10^^ 


0.486 (H 


-0.005 


-0.010) 


0.082 (H 


-0.020 


-0.015) 


0.78 


1.54 X IQi-* 


0.467 (H 


-0.012 


-0.014) 


0.073 (H 


-0.050 


-0.033) 


0.93 



the property that mass bins with a lower number of halos 
also h ave lower mean values, as in FigurelU lKasun fc Evrardl 
i2004f) find that the distribution about the mean "is well fit 
by a Gaussian" and contains no halos with low values of a. 
We find the same in our spherical window analysis. 

Now we turn our attention to the relationship between 
(s) and the dispersion for each bin. At z = we find a 
dependence of the dispersion on (s) (t op plot in Figure 13 . 
The relationship is steeper than that of ljing fc SutoH2002l) , 
who determined that as = 0.2 l(s) for the mass range they 
studied. However, one can see in the bottom plot of Figure 
13 that at z = 1 this relationship is no longer visible. This 
may be due to the fact that the number of halos in each 
bin at 2 = 1 is much lower and therefore dominated by 
systematics. It could also be that a large enough range in 
mass is not probed at z = 1 to see the relationship. We ar e 
unable to draw a conclusion similar to l.Tine fc Sutol l)2002l) . 
We therefore assume a constant value oi as =0.1, which is 
consistent with our results at all redshifts. 



4.5 Middle Axis Relationship 

The largest to smallest axial ratio, s, does not uniquely de- 
termine an ellipsoidal shape. There is still the determination 
of the relationship of the middle axis (6) to the s mallest (c) or 
larges t (a) axis. We find in our analysis, as did l.Tine fc Stitol 
l|2002), that the function P{p = c/h\s) exhibits a nice sym- 
metric behaviour. More commonly examined is the distribu- 
tion of g = b/a which can be trivially obtained from P{p\s). 
Figure |S| contains six histograms of p for different rang es in 
s. The curves are a fit proposed bv ljing fc Sutol i2002l) . 



P{p\s) = 



2(1 - ,5) 



1 - 



2p - 1 - 

1 - s 



with s - 



for S < Sn 



0.55 and s = s for s > Sn 



(9) 



It 



should also be noted that P{p\s) = below s. Ijing fc Sutol 
l|2002ll fit with a cut-off of Smin — 0.5, but otherwise we find 
agreement with their results. 



4.6 Radial Dependence of Shape 

The ellipsoidal shape of a halo is also found to be dependent 
on the radius at which the shape is determined. There is a 
systematic dependence of shape on radius with more massive 
halos having a steeper gradient in s with radius than lower 




0.6 0.65 



Figure 5. (s) versus CTs- Top: We find some evidence for a depen- 
dence of the dispersion on (s) at 2 = 0.0 if we perform a weighted 
linear least squares fit to the best values (black circles) of (s) 
a nd ais\ (black short dashed line) . It is slightly stepper than that 
of jjing fc Sulci ^2002^ ( green long dashed line). Also shown are 
the raw average and dispersion points (red triangles) Bottom: By 
2 = 1.0 this relationship seems to have disappeared. Due to the 
lack of a clear relationship between {s) and CTs, we favour a con- 
stant value of CTs = 0.1, wh ich is consistent wit h all redshifts and 
is roughly consistent with ljing fc Sutel i2002l) for the values of 
(s) probed. 
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Figure 6. The distribution of p = c/b in given bins of s show s 
a very similar behaviour to that found in | Jing fc Sutol J2OO2II 
The fi t line is from Equation originally found irt ljing^^Sutq 
J200I). 
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mass halos. In order to study the radial dependence of shape, 
halos in the L8O0.96 simulation were examined at 5 different 
fractions of their virial radius (Figure 0. For all halo mass 
bins there is a tendency for halos to be more spherical at 
larger radii, with more massive halos having a steeper change 
in (s) with radius. 

We also examined the value of (p) with radius and found 
no radial dependence. Therefore, a and b have the same ra- 
dial dependence and the largest axis a becomes relatively 
shorter with radius. We examine the relationship of the ra- 
dial dependence of s with mass and combine it with Equa- 
tion Q to find a shape-radius relationship, 

{s)(A/vir,r) = 6(MviO(r/rvi. - 0.3) + 0.54(Mvi./M*)-''-°' 

(10) 

with 



6(Mvir) = 0.0371ogio(Mvir/A/.) -I- 0.062. 



(11) 



We also examined the shape at different radii with a 
spherical window. The shape of the halo did not change 
very much o n average with radius. Th is is consistent with 
the result of iBailin fc Steinmetd l|200^ (see Section |SJ. 



4.7 Triaxiality 

Often ellipsoids are described in terms of their triaxiality 
(prolate, oblate, or triaxial). One way of expressing the tri- 
axiality of an ellipsoid is by using the triaxiality parameter 
iFranx. Ilhngworth fc de Zeeuwl (Il99ll) : 



T ■ 



(12) 



r/r . 

/ vir 



An ellipsoid is considered oblate it <T < 1/3, triaxial with 
1/3 < T < 2/3, and prolate ff 2/3 < T < 1. In Figure |S1 we 
divide up the halos into the same mass bins as in Figure |7| 
and analyse the triaxiality at i? = 0.3i?vir and JJvir- We find 
that most halos are prolate in shape with very few oblate 
halos, even at Rvii- The deficit of halos with T very close to 
1 is not physical. Due to the iterative process we use to de- 
fine shapes, if any two of the axes become degenerate (same 
length) the process has trouble converging. In most cases it 
does converge but with a large systematic error. Some au- 
thors have suggested that halos become oblate at large radii. 
We find only a small trend to less prolateness at large radii, 
but no evidence of a shift to oblate. Figure |S] also shows that 
larger halos, mainly those above M« , are almost entirely pro- 
late. Because we expect halos with masses above M« to be 
undergoing a higher rate of merging than halos with masses 
below M*, and because it has been sho wn that this merg- 
ing happens along preferred directions ( Kncbc ct al. 200^; 
FZentne r ct al. .2005b; F altenbacher et a l. 2005), the prolate- 
ness is most likely due to merging. This is in support of the 
idea that halo merging is responsible for the distribution of 
shapes. The internal velocity of a halo is also related to the 
merger history and therefore one would expect a relation- 
ship between he velocity structure of halos and their shape. 
We examine this in the next subsection. 



Figure 7. (s) as a function of radius at 2 = 0. 



4.8 Alignment with Velocity and Angular 
Momentum 



We have shown that the shape of a halo is related to the mass 
and have seen some clues that the shape is due to merging. 
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MerKing is also relate d to the angular momentum of halos 
jVitvitska et al.ll2002^ and their velocity dispersion. In or- 
der for a coUisionless system such as a DM halo to sustain 
its shape after merging, there must be an internal pressure 
provided by the velocity dispersion. If this is the case one 
would expect the internal velocities to be correlated with the 
shape. In order to investigate this we examine the alignment 
of the angular momentum and the velocity anisotropy of the 
halo with the shape. The angular momentum used here is 
calculated using the same particles found in the final ellip- 
soidal volume from our iterative method for determining the 
shape, although the results do not have a large dependence 
on which subset of par ticles within the halo is used. We find , 
as was pointed out bv lWarren et al.l l(l993),|^rme3 JlQQTt) . 
and subsequently seen by others, that the angular momen- 
tum is highly correlated with the smallest axis of the halo. 
In Figure|5]we show the absolute cosine of the angle between 
the indicated axis and the angular momentum vector. If the 
orientations were random the plot would be of a straight 
line at a value of 0.5. A peak at | cos6\ = 1 means that the 
axes are most often aligned and a peak at | cos 9\ = Q means 
that the axes are most often perpendicular to the angular 
momentum. As one can see, the smallest axis is most of- 
ten aligned and the largest axis is most often at an angle 
of 7r/2 from the angular momentum. Although the angular 
momentum is aligned with the smallest axis as would be 
expected for an object which is rotationally supported, DM 
halos are found not to be rotationally supported. Therefore, 
the significance of this alignment points not to a cause and 
effect relationship but to a shared origin. It has been shown 
in previous studies that the angular momentum of halos is 



0.2 0,4 0.6 0.8 1 

cos{fl) 

Figure 9. Top: The probability distribution of the cosine of the 
angles between the largest, middle, and smallest axis and the 
angular momentum vector. If the angular momentum were ran- 
domly oriented the graph would be a flat line at a value of 0.5. 
Bottom: The probability distribution of the cosine of the angles 
between the shape axes and velocity anisotropy axes. The velocity 
anisotropy is highly correlated with the shape. 



large ly determined by the last major merger jVitvitska et alJ 
I2OO2!), and that, at least during very active periods, merging 



fboth minor and major) happens along 


preferred directions 


jKnebe et al.l2nn4l:IZentner et a1.l2nn5bl: 


iFaltenbacher et alJ 


|2005|). It would seem, based on this, that the shapes and 



orientations of DM halos, at least during active merging pe- 
riods, can be attributed to directional merging. 

In order to determine whether halos are relaxed and 
self supporting we examine the relation of the velocity 
anisotropy to the shapes of halos. The velocity anisotropy is 
defined in the same way as the inertia tensor used to measure 
the shape. 



(13) 



We do not use a weighted version of this and we do not iter- 
ate, because neither of these make much physical sense. We 
calculate the velocity anisotropy tensor again using the par- 
ticles found within the ellipsoidal shell defined by the shape 
analysis. As with the angular momentum the alignment is 
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we determine the MAH for each halo at z = by followinR 
the ev olution of its most massive progenitors. IWechsler et all 
I2O.QJ) showed that the MAH of a halo can usually be well 
fit by a single parameter model, 



Figure 10. Triaxiality of the velocity anisotropy tensor and the 
shape for R = 0.3-Rvir (blue) and -Rvir (red). 



very insensitive to the particles used. We then determine the 
angle between the respective axis (i.e. a, b, & c). In Figure]^ 
we plot the distribution of absolute cosines between ashape 
and a„ei, bsUaps and &„ei, and Cshape and c^ei- From this one 
can see that all three of the axes are highly correlated. The 
strength of the alignment between the velocity anisotropy 
tensor and the shape suggests that the shape is supported 
by internal velocities. But how does the shape relate more 
directly to the velocity anisotropy? 

In Figure IIUI we show the triaxiality of the velocity 
anisotropy and the density. The velocity anisotropy on aver- 
age is more spherical in shape than the density. This is the 
expected trend from the Jeans equation for an ellipsoidal dis- 
tribution. The velocity anisotropy is directly related to the 
potential which has the same orientation as the shape but 
is more spherical due to the fact that potential is related to 
the spacial derivative of the density. It would therefore seem 
that the mass dependence of shape can not be explained 
by different relaxation times. This is also supported by the 
fact that halos are more aspherical in the centre where the 
relaxation time would be shorter than at the virial radius. 



5 MERGER HISTORY AND SHAPES 

So far we have investigated the evolution of halo shapes in 
fixed mass bins as a function of redshift and for two differ- 
ent values of as- Additional insight into the origin of shapes 
and their dispersion can be gained by tracing the evolution 
of individual halos. In order to quantify the evolution or 
mass accretion history (MAH) of the halos we have con- 
structed merger trees. Fo r more information on the merger 
trees, please see lAlkoodl (|200^. From these merger trees 



M(a) = Mo exp 



-2ac 



1 



- 1 



(14) 



where Mo is the mass of the halo a z = and ac is the scale 
fac tor at which the log sl ope of the MAH is 2. Although in 
the lWechsler et all i2002h they only allowed ac to be a free 
parameter, we find that by also allowing Mo to be a free 
parameter we are able to better recover ac for halos which 
had experienced a recent major merger. H alos with lower 
values of ac formed earlier, and as shown bv lWechsler et aD 
(2002), have a higher concentration. By means of Equation 
1141 1 we assign an ac to every halo found at z = 0. 

In Figure 1111 we plot s versus ac for the halos in the 
LI2O0.9 simulation split into separate mass bins. We find 
that halos which formed earlier are on average more spheri- 
cal with a dispersion of 0.08 — 0.1 (see Figure lT^ for all mass 
bins. This implies that the scatter in the {s)(Mvir) relation 
can not be completely attributed to the different values of 
ac for that particular mass bin. However, the dependence on 
ttc is less for the higher mass halos and this would explain 
the mean-dispersion relationship explored in Subsection l4.4l 
Higher mass halos were found to have a smaller dispersion 
than lower mass halos at z = 0. This can also be seen in 
Figure m\ It is very likely that the residual scatter is due, 
at least in part, to the pattern of infall. Since s is derived 
from an inherently three dimensional quantity, namely the 
inertial tensor, a one dimension parameterisation may not 
be sufficient to capture all of the physics involved. A more 
careful study of infall is needed to explain the dispersion in 
s. 

The above investigation makes clear that the dispersion 
of halo shapes cannot be explained by appealing to a single 
parameter description of the MAHs. However, an average 
evolutionary pattern for halos which is dependent on both 
ac and mass is seen. Figure 1131 displays the evolution of 
(s) (sorted by ac) with scale factor a for a particular mass 
bin at z — 0. Halos that formed early (lower ac) are more 
spherical today as was pointed out above. Moreover, they 
become spherical more rapidly (indicated by the increasing 
slopes for halos of low ac), although the transformation rate 
towards spherical shapes seems to slow for all values of ac 
with increasing expansion factor. In Figure 1131 the results 
for the lowest mass bin (3.2 x 10^^ < Mvir < 6.4 x 10^^) are 
shown. Apart from a systematic shift to lower values of s the 
corresponding plots for higher mass bins look very similar. 

Halos which have early formation times (low ac) at a 
fixed mass today have typically accreted more mass since ac 
than halos with higher values of ac. The rapid transforma- 
tion towards spherical shapes for early forming halos implies 
either that lower mass halos become spherical more rapidly 
after ac, probably due to shorter dynamical times, or that 
mass accretion after ac is more spherical, therefore causing 
the halo to become more spherical as well. By examining 
other mass bins we find that halos of different masses today 
but with comparable values for ac, thus different masses at 
ac, show the same rate of change in s, but with different ini- 
tial values of s. This finding suggests that the rate at which 
a halo becomes spherical depends on its Oc rather than on 
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Figure 11. (s) vs characteristic formation epoch for different 
mass bins (mass quoted in units of 10^^ Mq). Only bins that 
contain at least 10 halos are shown (square points). There is a 
distinct trend of shape with ac for the lower mass bins. At higher 
mass there is still a trend but it is uncertain how strong the trend 
is due to the lower number statistics. Solid black line is a linear fit 
to the points and dashed line is the la scatter about the points. 

its mass. We find that we can approximate the dependence 
of s on the expansion factor a for a > Qc 0.1 by a simple 
power law 

(s) (a) oc (a — Oc)" : a >ac + 0.1, (15) 

where v has to be fitted for the particular halo. In Figure 
1141 we display the values of v versus flc determined by 
fitting for the LSOo.gb simulation. The L8O0.96 was divided 
into bins of log mass and ac- The average MAH for each 
bin was fit by Equation 11511 using minimisation. All bins 
containing at least 20 halos were used determine the function 
V. We find a tight correlation between the u and ac which 
can be approximated by 

u = 1.74 X aZ°-^ . (16) 

This fit is represented by the solid line in Figure 1141 The 
remarkable success of Equation 11511 to fit the data supports 
the idea, that the transformation from aspherical to spheri- 
cal halo shapes is driven by mass accretion becoming more 
spherical after ac- The physical reason for the observed be- 
haviour merits further investigation. 



6 COMPARISON WITH PREVIOUS 

DETERMINATIONS OF HALO SHAPE 

In the previous sections we have explored many aspects of 
halo shapes. Central to this discussion has been {s)(Afvir). 



Figure 12. Distributions of s for a given range in ac and mass. 
The solid histogram is for the mass range: 3.2 X 10^^ < Mvir < 
6.4 X 10^^ and the dotted histogram is for a mass range of 1.28 X 
10^3 < Mvir < 5.12 X 10^3. The dispersion m s m a given mass 
bin can be explained in part by the different MAHs. 

This relationship has been examined by many recent studies 
all of which seem to determine different relationships. In this 
section we address these discrepancies. 

First, an examination of the difference in the inferred 
shape from the use of the weighted versus the unweighted 
inertia tensor is needed. Most recent authors prefer the 
weighted (or reduced) inertia tensor (Equation 0) which 
is the method we have chosen to use. The motivation for 
the use of the weighted inertia tensor, /, is due to the bias 
present in the unweighted method to particles at larger radii. 
By weighting the contribution from each particle in the sum 
by the distance to the particle squared, / is less sensitive to 
large substructure in the outer regions of the analysis vol- 
ume. To test the difference between the methods we exam- 
ined a sample of halos from the LI2O0.9 simulation using the 
iterative method with both versions of inertia tensor (Figure 
I15II . Both iterative methods give similar results for the mean 
quantities (inset in Figure [T^ as a function of mass with in- 
dividual halos differing by As < 0.15. The detailed distribu- 
tions are different and have some interesting features (Table 
Ol. We find that halos which have lower s values for the 
weighted method over the unweighted method have larger 
substructure near the centre and halos which have lower s 
values for the unweighted method had larger substructure 
near the outer edge of the analysis volume. At very small 
axial ratios the unweighted method seems to always give 
larger axial ratios. As we have shown, halos with s < 0.3 are 
late forming and are strongly contaminated by substructure. 
This leads to the unweighted method giving a lower value 
of s indicated by the high skewness shown for the method 



12 Allgood et al. 




Figure 13. The evolution of halos with different values of Oc in 
the mass bin 3.2x 10^^ < M^i^(z = 0) < 6.4x 10^^. Halos become 
more spherical after a short period after ■ The halos which form 
earlier become spherical more rapidly. Log binning was chosen to 
even out the number of halos in each bin. 



Figure 14. The rate of change exponent v (see Equation 1151 '! 
versus expansion factor at halo formation ac- The solid line dis- 
plays the fitting formula given by Equation 1161 . Only bins con- 
taining at lease 20 halos is displayed. The ac value of each point 
is the average value for that respective bin and the errorbars rep- 
resent the variance determined by the fitting. 



Table 3. Weighted vs. Unweighted Shapes 



Method Mean Mass (s) as kurtosis skewness 



weighted 


5 


X 


10" 


0.694 


0.094 


-0.029 


-0.132 


unweighted 


5 


X 


10" 


0.686 


0.097 


-0.302 


-0.053 


weighted 


5 


X 


10" 


0.614 


0.111 


0.253 


-0.188 


unweighted 


5 


X 


10" 


0.603 


0.118 


0.605 


-0.287 


weighted 


5 


X 


10" 


0.518 


0.117 


1.756 


-0.532 


unweighted 


5 


X 


10" 


0.524 


0.126 


1.628 


-0.592 


weighted 


5 


X 


10" 


0.455 


0.114 


2.066 


-0.515 


unweighted 


5 


X 


10" 


0.472 


0.133 


1.428 


0.150 



in Table El by the skewness. For the two well resolved mass 
ranges (10^^ < Mvir < lO" and lO" < Afvir < lO") the un- 
weighted method has a larger negative skewness. Note that 
all mass bins except for the last unweighted bin have nega- 
tive skewness (this was discussed in Section|3J, but the most 
massive bin suffers from low statistics. The distribution of 
as values is always broader for the unweighted method. 

There are more differences than just the form of the 
inertia tensor used. In order to compare our results to a 
selected number of previous results (Figure [T^ we have re- 
peated the shape analysis using the methods described in 
the corresponding papers. The previo us work which i s mos t 
similar to the current work is that of ISoringel et all i2004l) . 



Our findings are very similar to the ISpringel et al.l (|200^ 

results except we found that (s)(Mvir) has a slightly higher 
normalisation. Through private communication with Volker 
Springel, we were provided with an updated set of data 
points which come from a more complete sample and are in 
much better agreement with our results (open green squares 
in Figure Not only do our results agree at z = 0, but 
also at higher redshift (see Figure 

IJing fc Sutcl ll2002l) ( JS) studied 12 high resolution clus- 
ters with A'^ ~ 10® particles and five cosmological sim- 
ulations with A'' — 512^ particles in a lQQh~^ Mpc box 
with both an SCDM and ACDM cosmology. The simula- 
tions were performed with a P^M code with fixed timestep- 
ping and a spatial resolution of 10 — 20/i~^ kpc. They used 
a FOF halo finder and analysed the shapes of the high 
resolution clusters in isodensity shells as a function of ra- 
dius, finding that the halos are more spherical at larger 
radii. After determining the relationship of shape with ra- 
dius they developed a generalised ellipsoidal NFW den- 
sity profile (Figure 1171 shows that our method gives sim- 
ilar results). They applied this generalised fitting routine 
to the cosmological simulations and determined generalised 
NFW parameters and shapes statistics. The shapes were 
determined using an isodensity shell at an over density of 
2500pc (where is the critical density) which corresponds 
roughly to J? = 0.37?vir. The mass range analysed only cov- 
ered one order of magnitude in mass (2.1 X W^h-^ Mq < 
Mvir < 1 X 10"/i~^ Mq). They found a result very simi- 
lar to ours, {s)(Mvir) = 0.54(Mvir/M.)-"-''''["'"(^)l°'\ and 
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Figure 15. Axial ratios for lialos in a cosmological simulation 
are divided into mass bins and tlie shapes are calculated using 
the weighted (s^,) and non-weighted iterative inertia tensor 
methods. The two methods agree within ^ 10% and give the same 
value when averaged over a given mass bin (inset graph). 

dispersion which is well fit by a Gaussian distribution with 
CTs(Mvir) = 0.113(Mvi^/M.)-°•°''['^'"(^)l'' ^ We do not find 
any evidence for a steepening of the exponent with redshift 
as they do. 

The disagreement between our findings and those of JS 
regarding the scaling of (s) with mass is due to the proce- 
dures used. In the JS analysis they determine the shape 
of an isodensity shell at 2500pc ± 3%, completely ignor- 
ing the interior of the shell. JS analysed halos with masses 
greater than 6.2 x lQ^^h~^ M© which tend to be dynam- 
ically young and often have double cores. This can affect 
the shape a lot, but their analysis would not pick this up, 
due to the neglect of the central region. In our iterative 
inertia tensor analysis we include the centres. In order to 
confirm that the difference is truly due to the shell versus 
the solid ellipsoid we analysed halos from the L200o.g sim- 
ulation in the mass range 1 — 4 x 10^^/i~^ using the 
technique presented in JS. We examined isodensity shells 
at 2500pc with a thickness of ±30%, instead of the ±3% 
used by JS. We needed to examine thicker shells in order 
to obtain enough particles to do the analysis because the 
L2OO0.9 has less mass resolution than the simulations stud- 
ied by JS. We found that the inertia tensor method gives 
{s)tensor = 0.485 ± 0.008 aud as = 0.091 ± 0.006 and the JS 
method gives (s)js = 0.515 ± 0.008 with the same scatter 
for the same halos. The difference is due to the fact that a,js 
is systematically larger at low stmsor- This pattern is born 
out by a quantitative analysis. When we split the sample 
at Stensor = 0.45 (roughly where agreement begins) , we get 
that above Stensor ~ 0.45 the samples agree quite well with 



Figure 16. Comparison of (s)(Mvir) relation with previous stud- 
ies. We attempt to reconcile the differences between our results 
and those of other authors. We present the results of a shape anal- 
ysis of the L80o.9i,, L120o.90r, and L200o.9 simulations using the 
iterative inertia methods at R = 0.3-Rvir (black, pink, cyan) and 
non-iterative spherical window analysis at i?vir (red, green, violet). 
In addition, we present the results of a shape analysis of the 
L80o.9i) and L120o.90r using the iterative method at -Rvir (blue). 
The black line is our proposed fit from Equation and this 
should be compared to the resul ts of Springel (priva te communi- 
cation) (green open squares) and l.Tinef fc Sutol i2n02ll (orange dot 
dash). The blue line is a fit to the blue points, which should be 
compared to the lBuUockl l2002h line (violet long dash). Finally the 
red line is a renormalised version of the|Kasun & Evrard |20o3) 
fit which should be compared their fit (brown small dash). The 
thin black dot dash line a t (s) ~ 0.7 is the spherical shell fit of 
iBailin fc SteinmetJ l200.'t) . The bold portions of the lines indicate 
the mass ranges where the fit was compared to simulated data by 
the respective authors. 



{s)tensor = 0.571 ± 0.007 vs (s) JS = 0.576 ± 0.011. Whereas 
below Stensor = 0.45 wc get {s)tensor = 0.424 ± 0.007 vs 
{s)js = 0.472 ± 0.008. The diflterence at the low end is due 
to the missing of the dynamically active cores by JS. If JS 
had extend their analysis to lower mass halos were multiple 
cores are not as common their determination of (s) would 
converge with ours. Because X-ray observers normally do 
not choose to only analyse the outer shells of clusters due 
to the fact that the X-ray observations get noisier with the 
distance from the centre and because optical observers may 
not see the multiple cores when analysing cluster member 
velocities, we prefer the method which includes the effect 
of the multiple peaks. In Paper II we show that using our 
method with some additional assumptions one can account 
for the observed X -ray ellipticity measurements. 

iBulloclj (|20o3) analysed the shapes of halos in a ACDM 
simulation with erg = 1.0 at three different redshifts {z = 
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Figure 17. Comparison of an isodensity shell (red) and a tensor 
ellipsoid (blue). Particles are selected by the JS (red) and inertia 
tensor (blue) methods and projected on the x-y plane of the sim- 
ulation box. The shortest/longest axis ratio is s = 0.49 (0.48) for 
the isodensity shell (tensor ellipsoid). The semi-major axis of the 
isodensity shell 0.23Rviri consistent with JS for 2500pc isodensity. 



agreement. In examining the other two simulations for lower 
mass halos we are unab le to recover the extrapolation of the 
iKasun fc Evrar d* f2004*) relationship. In fact, we see a tran- 
sition from the Kasun_ fc Evrard ( 2004) relationship to the 
relationship of Ba ilin fc Steinmet jH2005h (discussed below). 
We also find that the mean shape relationship has almost no 
dependence on radius when using a spherical window func- 
tion . 

iBailin fc Steinmetd (jioO^ analyse the shapes of halos 
at different radii in spherical shells in the mass range of 
10"/i"^ M0 < M < 5 X 10". After determining the axial 
ratios they then apply an empirical correction of Strue = 
•Smeasure to correct for the use of a spherical window. They 
find that all halos have an axial ratio of (s) ~ 0.63 at 
7? = 0.4_Rvir with the scaling applied, which implies that 
they measure (s) ~ 0.766 in their spherical window. This 
result is in very good agreement with our spherical analysis 
(green and red data points in Figure lTSl . However, we do not 
find that halos of different masses have the same mean axial 
ratio. There seems to be some evidence that the {s)(Mvir) 
relationship fiattens out below M», but it is definitely not 
constant with radius. Simulations with even higher mass res- 
olution are needed to investigate for the possibility of fiat- 
tening below M«. For an extra check we also analysed the 
halos in a spherical shell between 0.25i?vir and 0.4i?vir and 
measure a roughly fiat value for all halos of s = 0.77. The 
disagreement about the (s)(Af) relationship most likely lies 
in the determination of the empirical spherical window cor- 
rection. The correction was determined using Monte Carlo 
halos with no substructure, but we find that substructure 
plays a role in the determined shape of the halo. 



0.0, 1.0, 3.0). The simulations were performed using the ART 
code in a 60/i~^ Mpc box with 256^ particles and spatial 
resolution of \Mi~^ kpc. The analysis of shape was done 
using the weighted inertia tensor in a spherical window with 
K = -Rvir. The axial ratios were determined iteratively un- 
til convergence was obtained using a similar criterion as we 
have used, but the window remained spherical. The use of 
the weighted inertia tensor and iterative axial ratio deter- 
mination seemed to almost eliminat e the effect of using a 
spherical window (discussed below). iBulloclj ll2002l) found 
that (s)(Mvir) =i 0.7(Mvir/10i2/i-i Mp))-°-°'^(l-Fz) -°-^ was 
a good fit to the simulation. The empirical scaling of lBulloc^ 
i2002|) is similar to what we find, but the powerlaw is 
steeper. This can be attributed to the lower resolution and 
possibly the use of a spherical window. Bullock's higher nor- 
mal isation is due to the hig her as- 

iKasun fc EvrardI ll2004l) determined the shapes of clus- 
ter halos (M200 > 3 X VS'^hr^ Mq) in the Hubble Vol- 
ume simulation. They calculated the axial ratios using the 
unweighted inertia tensor in a spherical window at -R200, 
the radius of the sphere within which the mean density 
is 200pc(-z), with pc{z) being the critical density at red- 
shift z. They determined a relationship of (s)(Mvir) = 
0.631[1- 0.023 ln(Mvir/10"/i"^ M0)](l-H2)~° °**^ We com- 
pare our analysis with theirs by performing the same spher- 
ical analysis at i? = -Rvir, which is slightly larger than -R200- 
We find that in our largest box simulation (L2OO0.9) where 
we have good statistics on cluster mass halos we find good 



7 COMPARISON WITH OBSERVATIONS 

Since all of the differences between the shape statistics ex- 
tracted from pure collisionless simulations by various au- 
thors can be reconciled by considering the different methods 
used to determine shapes, a comparison between observa- 
tions and simulations is in order. Much of the attention halo 
shapes have received lately is due to the recent estimates 
of the shape of the Milky Way's host halo. Most estimates 
find the Milky Way's host halo to have an oblate shape with 
s > 0.8. This is in contrast with s 7^ 0.6±0.1 for lO^^/i"^ M© 
halos found in pure collisionless simulations, though there 
is some evidence that the halos become more spherical 
when ba ryonic cooling is included iKazantzidis et al]|2004l : 
iBailin et al„, .2Q05i) and that some become oblate. The pres- 
ence of gas cooling will invariably make the halos more 
spherical but the extent of the effect is not yet fully under- 
stood. Recently, there have however been studies of the M 
giants in the leading edge of the Sagittarius dwarf stream 
t Helrni 2004; Law ct al. 2005), which concluded that the 
best fit shape of the host halo is a prolate ellipsoid with 
s = 0.6. 

Another way of measuri ng the shape o f DM ha- 
los is through weak lousing. iHoekstra et alJ (|2q3) S'ld 
iMandclbaum et al. (2005) performed studies of galaxy- 
galaxy weak lensing using the Red-Sequence Cluster 
Survey and the Slo an Digital Sky Survey respectively. 
iHoekstra et al.l (|22q3) determine the average shapes of halos 
by measuring the orientation of the galaxies, then stacking 
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Figure 18. Projected ellipticity of galaxy mass halos at 2 = 0.33 
for (Tg = 0.9 (black solid) and erg = 0.75 (red dashed). 

the galaxy images with the orientations ahgned, and finally 
measuring the shear field around the stacked image. This 
measurement of the shear provides a rough estimate of the 
dark matter halo shapes at z ~ 0.33. They found an av- 
erage projected ellipticity of (e) = (1 — q2D) = 0.20tg;Q5, 
corresponding to s = 0.66ln' nL for halos with an average 
mass of 8 X 10"/i-i Mq. The lHoekstra et alJ (|200i) deter- 
mination was hindered by the fact that the galaxies were 
stacked together regardless of morphological type, and one 
would reasonably assume that morphological type is related 
to merger history and possibly orientation with the host 
halo which will in turn affect the measured halo shapes as 

we have sho wn. 

1-1 the iMandelbaum et alJ l^20^!^^ analysis of SDSS 
galaxies colour was used a s a proxy for morphological type. 
IMandelbaum et al.l l^200!^^ studied 2 million lens galaxies 
with r > 19 and 31 million source galaxies dividing the 
lenses into bins of colour and luminosity. They find a sug- 
gestion that spiral (blue galaxy) light ellipticities may be 
anti-aligned with the halo ellipticities at the 2a level and 
a suggestion that elliptical (red galaxy) host halo elliptic- 
ities are luminosity dependent. Since we and others have 
shown that halo angular momentum is highly aligned with 
the smallest axis, this finding would suggest that the an- 
gular momentum of spiral galaxies is not aligned with the 
angular mo mentum of th e their host halos. Recent theoret- 
ical work bv lBailin et alJ (|200^ seems to support this idea 
as well, although they only study eight spiral galaxy simu- 
lations. If this were t he case one would assume that in the 
iHoekstra et alJ (|200^ analysis the measured signal would be 
diminishe d by this. It would no t be completely nullified be- 
cause the lHoekstra et all l|2(^(j 4[l sample is domina ted by el- 
liptical galaxies. Indeed. IMandelbaum et al.l (120051) show by 



combining the appropriate luminosity and colour bins that 
their fi ndings are in agreement with those of iHoekstra et all 
('2004") . In Figure we have plotted the distribution of pro- 
jected axial ratios for the L80o.9a and L8O0.75 simulations 
for 1000 random lines of sight through the box for galaxy 
mass halos. In projection the differences between the simula- 
tions become small hindering any sort of determination of erg 
via lensing studies. The peaks of the distributions a re both 
broad ly in agreement with th e findi ngs of .Hoekstra et alJ 
J2004I) and IMandelbaum et all feOOSl) . Further study of the 
galaxy / host halo alignment relationship and galaxy mor- 
phology / halo merger history relationship is in order to bet- 
ter understand and predict the galaxy-galaxy lensing mea- 
surements. 



8 SUMMARY AND DISCUSSION 

First we investigated the dependence of halo shapes on their 
masses. Our halo sample, based on six different simulations, 
covers three orders of magnitude in mass from galaxy to 
cluster scales. Furthermore we have examined halo shapes 
as a function of redshift and ag. Our analysis of shape is 
based on the halo volume enclosed hy R ~ 0.37?vir. This 
particular radius is chosen for several reasons; halos should 
be fairly well relaxed within this radius and shapes measured 
from X-ray gas in cluster should extend to this radius. At 
this and other radii we find that the mean shape of dark 
matter halos depends on the halos mass. We find that the 
mean largest-to-smallest axial ratio s — c/a at radius 0.3i?vir 
is well described by 

{s){AU,,, z) = (0.54 ± 0.03) (17) 

The distribution of s in each mass bin is well fit by a Gaus- 
sian wi th a = 0.1. The rel ation found here is steeper than 
that of Ijing &: Sutd i2002h at 2; = 0, thus predicting less 
spherical cluster mass halos and more spherical galaxy mass 
halos. 

In order to reconcile our results with sometimes strongly 
deviating findings by other authors we have applied their 
particular methods. We find that the disagreements be- 
tween different studies are mostly due to differences be- 
tween methods of measuring axial ratios. However, there 
remains one open question brought up by this comparison. 
T here seems to be an inc onsistency between the analysis 
bylKasun & Evrard^ ('2004'), if extended to galaxy mass ha- 
los, and Bailin & Stcinmctz (200^. Both groups analysed 
halo shapes within a spherical window. We find that the ap- 
plication of a spherical window leads to a transition from 
a power-law like behaviour above h'h to a mass indepen- 
dent shape below M, . This transition is found in the gap 
between the mass ranges analysed by the two groups and 
therefore was not recognised by either of them. Possibly a 
similar transition can be found applying our method of an it- 
erative ellipsoidal window. But the effect seems to be smaller 
and requires, if apparent, even higher mass resolution than 
utilised in this work. 

We find that the mean shape relation becomes shal- 
lower and more spherical at increased radius. We also find 
that higher mass halos have a steeper relationship between 
shape and radius than smaller mass halos. Since cluster sized 
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halos are on average younger than galaxy sized halos, we are 
comparing dynamically different objects. The presence of an 
increased amount of massive substructure near the centre of 
dynamically young objects may be the reason for the steeper 
relation of shape with radius for cluster mass halos than 
galaxy mass halos. 

Our analysis of the halo shapes as a function of redshift 
leads to the following results. Within fixed mass bins the 
redshift dependence of {s)(Mvir) is well characterised by the 
evolution of M*, unhke the findings of I.Ting fc Sutol ()2002f ) 
who predict a much steeper relation of oc M^^"^ at high 
redshifts. We find that Equation H17^ works well for different 
values of as (a variation of results in a variation of M» 
which appears as a normalisation parameter in the {s)(Mvir) 
relation). Also worth noting is that at higher redshift the 
possible broken power-law behaviour disappears, but if it 
were truly due to M* we would expect this, because already 
hj z = 1 Mf is below our mass resolution. 

We find that the {s)(Mvir) at z = for galaxy mass 
halos is ~ 0.6 with a dispersion of 0.1. This result is in 
good agreement with only one estimate for the axial ra- 
tio of the Milky Way (MW) halo. Helmi (.2004} claims 
that a study of the M giants in the leading edge of the 
stream tidally stripped from the Sagittarius dwarf galaxy 
leads to a best fit prolate halo with s = 0.6. After anal- 
ysis of the same d ata. iLaw et al . (2005) confirm this find- 
ing. Other studi es ||lbM^^ ^ni2001 : Maiewski et al. 2Di}3; 
iMartfnez-Delgado et al.ll2004h which examined different as- 
pects of these streams concluded that the MW halo is oblate 
and nearly spherical with q > 0.8. If the shape of the MW 
halo is truly this spherical, it is either at least 2a more spher- 
ical than the median, or else baryonic cooling has had a 
large effect on the shap e of the dark matter halo (see e.g., 
iKazantzidis et al.ll20oi) . 

Describin g halo shapes by t he triaxiality parameter T 
introduced bv lFranx et al.l (^^), we find that the majority 
of halos are prolate with the fraction of halos being pro- 
late increasing for halos with Mvir > M,. Since halo shapes 
are closely connected to their internal velocity structure, we 
compute the angular momentum and the velocity anisotropy 
tensor and relate them to both the orientation of the halo 
and the triaxiality. In agreement with previous studies we 
find that the angular momentum is highly correlated with 
the smallest axis of the halo and that the principal axes of 
the velocity anisotropy tensor tend to be highly aligned with 
the principal axes of the halo. The strong alignment of all 
three axes of the two tensors is remarkable since the velocity 
tensor tends to be more spherical, thus the determination of 
its axes might be degenerate which would disturb the cor- 
relation with the spatial axes. If the accretion of matter 
determines the velocity tensor the tight correlation between 
velocity tensor and density shape argues for a determination 
of the halo shape by merging processes. 

Finally we examine the evolution of shapes by following 
the merger trees of the individual halos. We find that the dif- 
ferent mass accretion histories of halos cannot fully explain 
the observed dispersion about the mean s within fixed mass 
bins. It is likely that an analysis of the three dimensional 
accretion is essential for the explanation of the dispersion 
at a fixed value of mass and a^. However, halos with ear- 
lier formation times (lower Uc) tend to be more spherical at 
z = 0. Furthermore, there is a pattern of halos becoming 



spherical at a more rapid rate for halos that formed earlier 
and this rate appears to be independent of the final mass. 
The evolution of the shape for a > Qc -1-0.1 is well described 
by 

(s)(a) oc (a-flc)" , (18) 

where u — 1.74 x a^"''^. We detect a definite trend for the 
transformation from highly aspherical to more spherical halo 
shapes after Uc- The change of s seems to be less dependent 
on the total halo mass but strongly infiuenced by the relative 
mass increase since which suggests that halos are becom- 
ing more spherical with time due to a change in the accretion 
pattern after from a directional to a more spherical mode. 
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APPENDIX A: RESOLUTION TESTS 

A potential source of systematic error in the determination 
of halo shapes is the limited number of particles involved for 
low mass halos. In Figure I^Tl we show the result of a Monte 
Carlo (MC) test. Our dark matter halos are adequately de- 
scribed by an elliptical NFW density profile independently 
of cosmological model. Therefore, in the figure we show the 
result of applying the reduced tensor method of Equation 
^ to MC halos built to have an ellipsoidal NFW profile. 
The axial ratios of the MC halos are drawn from Gaussian 
distributions of mean (s) = 0.7 and (q) — 0.85, and dis- 
persion a = 0.1. Approximately 450 halos were generated, 
each having ~ 1000 particles, in order to have a catalogue 
comparable to the sample of halos in the L8O0.9 box with 
mass in the range A/vir = lO"'^ - 10"-^/i"^ Mq. The scatter 
plot shows that individual values of s determined at 0.3i?vir 
which contains roughly 300 of the 1000 particles in the halo. 
The recovered shape can be in error by up to ~ 0.1. However, 
the scatter and mean of the distribution are very well de- 
termined by the inertia tensor method. The inset shows the 
histogram for the input values of the MC halos (solid line), 
and the histogram for the output values (i.e. determined by 
the tensor method; dotted line). Therefore, we conclude that 
the tensor method underestimates s by only 0.03 for halos 
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Figure Al. Results of applying our shape determination proce- 
dure at 0.3-Rvir to 450 Monte Carlo halos produced with deter- 
mined axial ratios. We found that the error in the recovered value 
of s could be as large as ~ 0.1. 
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Figure A2. (s) versus mass. This plot is a replica of Figure [Tl 
except we show mass bins below the determined resolution limit. 
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Figure A3. A direct comparison of the axial ratios of the most 
massive halos in the L102o.9r simulation to the corresponding 
halos in the LI2O0.9 simulation. 

of this mass. For halos of mass Mvir ~ 10"-^/i"^ Mq, the 
error falls to 0.01. 

The one problem with the above test is that we are 
attempting to recover shapes from smooth halos. In reality 
halos have substructure and the substructure must play a 
role in the shape of the halo. The presence of dense lumps 
close to the core will bias (s) to lower values relative to val- 
ues determined from isodensity shells. In order to test our 
ability to recover the axial ratios in a cosmological simu- 
lation we determine the mean shape of halos down to very 
low masses for all of our simulations and then compare them 
to one another. The main difference between this test and 
the previous one is that the halos in the cosmological sim- 
ulations contain substructure, but we do not know a priori 
what the distribution of shapes should be. In Figure IX2I we 
have plotted the shapes of halos down to very low parti- 
cle numbers. There appear to be two resolution effects at 
work here. The first effect is an extension of the result we 
found in the MC test above. At small particle number the 
recovered shape becomes very aspherical and all of the sim- 
ulations turn over at the same number of particles. For all 
of the simulations the turnover is detected at ~ 3000 par- 
ticles within _Rvir (which is ~ 1000 particles within 0.3-Rvir, 
where we are determining the shape). This is consistent with 
our MC tests. At higher particle numbers there seems to be 
another effect driven by the particle number and not the 
mass. Halos containing rip < 7000 particles for any given 
simulation show a trend of becoming more spherical on av- 
erage than the simulations of higher resolution for the same 
mass. In Figure IX3I we show the determined value for s for 
the most massive halos shared by the LI2O0.9 simulation 
and the resimulated subregion, L120o.9r. We find that halos 
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shapes between the two simulations can differ by as much as 
0.15 in s. This is not unexpected because in the resimulation 
process the identical halos are captured at slightly different 
times. The main point here is that the average value of s is 
lower in the higher resolution simulation for the same halos 
by ~ 0.05. 



REFERENCES 

Allgood B., 2005, PhD Thesis: Shapes and Orientations of 

Dark Matter Halos 
Bailin J., Kawata D., Gibson B., Steinmetz M., Navarro J., 
Brook C, Gill S., Ibata R., Knebe A., Lewis G., Okamoto 
T., 2005, ApJL, 627, L17 
Bailin J., Steinmetz M., 2005, ApJ, 627, 647 
Barnes J., Efstathiou G., 1987, ApJ, 319, 575 
Bryan G. L., Norman M. L., 1998, ApJ, 495, 80 
Bullock J. S., 2002, in Proceedings of the Yale Cosmology 
Workshop "The Shapes of Galaxies and Their Dark Mat- 
ter Halos", New Haven, Connecticut, USA, 28-30 May 
2001.Priyamvada Natarajan, Ed.. Singapore: World Sci- 
entific, 2002, |astro-ph/0106380| Shapes of dark matter 
halos. p. 109 
Buote D. A., Xu G., 1997, MNRAS, 284, 439 
Dubinski J., 1994, ApJ, 431, 617 
Dubinski J., Carlberg R. G., 1991, ApJ, 378, 496 
Faltenbacher A., Allgood B., Gottloeber S. amd Yepes G., 
Hoffman Y., 2005, MNRAS in press fastro-ph/0501452") 
Flores R. A., Allgood B., Kravtsov A. V., Primack J. R., 
Buote D. A., Bullock J. S., 2005, MNRAS submitted 
astro-ph/0508226 
Franx M., lUingworth G., de Zeeuw T., 1991, ApJ, 383, 112 
Frenk C. S., White S. D. M., Davis M., Efstathiou G., 1988, 

ApJ, 327, 507 
Helmi A., 2004, ApJL, 610, L97 

Hoekstra H., Yee H. K. C, Gladders M. D., 2004, ApJ, 
606, 67 

Hopkins P. F., BahcaU N. A., Bode P., 2005, ApJ, 618, 1 
Ibata R., Lewis G. F., Irwin M., Totten E., Quinn T., 2001, 
ApJ, 551, 294 

Jing Y. P., Mo H. J., Borner G., Fang L. Z., 1995, MNRAS, 
276, 417 

Jing Y. P., Suto Y., 2002, ApJ, 574, 538 

K asun S. F., Evrard A. E., 2004, ApJ submitted 

( |astr o-ph/0408056 ) 
Katz N., 1991, ApJ, 368, 325 

Kazantzidis S., Kravtsov A. V., Zentner A. R., Allgood B., 

Nagai D., Moore B., 2004, ApJL, 611, L73 
Klypin A., Kravtsov A. V., Bullock J. S., Primack J. R., 

2001, ApJ, 554, 903 
Klypin A. A., Gottlober S., Kravtsov A. V., Khokhlov 

A. M., 1999, ApJ, 516, 530 
Knebe A., GiU S. P. D., Gibson B. K., Lewis G. F., Ibata 

R. A., Dopita M. A., 2004, ApJ, 603, 7 
Kolokotronis V., Basilakos S., Plionis M., Georgantopoulos 

I., 2001, MNRAS, 320, 49 
Kravtsov A. V., Berlind A. A., Wechsler R. H., Klypin 

A. A., Gottlober S., Allgood B., Primack J. R., 2004, ApJ, 

666, 

Kravtsov A. V., Klypin A. A., Khokhlov A. M., 1997, ApJS, 
111, 73 



Law D. R., Johnston K. V., Majewski S. R., 2005, ApJ, 
619, 807 

Lee J., Suto Y., 2003, ApJ, 585, 151 

Majewski S. R., Skrutskie M. F., Weinberg M. D., Os- 

theimer J. C, 2003, ApJ, 599, 1082 
Mandelbaum R., Hirata C. M., Broderick T., Sel- 

jak U., Brinkmann J., 2005, MNRAS submitted 

Castro-ph/0507108 ) 
Martmez-Delgado D., Gomez-Flechoso M. A., Aparicio A., 

Carrera R., 2004, ApJ, 601, 242 
McMillan S. L., Kowalski M. P., Ulmer M. P., 1989, ApJS, 

70, 723 

Merrifield M. R., 2004, in lAU Symposium The Galactic 

Halo and CDM. p. 431 
Mohr J. J., Evrard A. E., Fabricant D. G., Geller M. J., 

1995, ApJ, 447, 8 
Navarro J. F., Frenk C. S., White S. D. M., 1996, ApJ, 462, 

563 

Oiling R. P., Merrifield M. R., 2000, MNRAS, 311, 361 
Springel V., White S. D. M., Hernquist L., 2004, in lAU 

Symposium The shapes of simulated dark matter halos. 

p. 421 

Tasitsiomi A., Kravtsov A. V., Gottlober S., Klypin A. A., 

2004, ApJ, 607, 125 

Thomas P. A., Colberg J. M., Couchman H. M. P., Efs- 
tathiou G. P., Frenk C. S., Jenkins A. R., Nelson A. H., 
Hutchings R. M., Peacock J. A., Pearce F. R., White 
S. D. M., 1998, MNRAS, 296, 1061 

Tormen G., 1997, MNRAS, 290, 411 

Vitvitska M., Klypin A. A., Kravtsov A. V., Wechsler 

R. H., Primack J. R., Bullock J. S., 2002, ApJ, 581, 799 
Warren M. S., Quinn P. J., Salmon J. K., Zurek W. H., 

1992, ApJ, 399, 405 
Wechsler R. H., Bullock J. S., Primack J. R., Kravtsov 

A. v., Dekel A., 2002, ApJ, 568, 52 
Zentner A. R., Berlind A. A., Bullock J. S., Kravtsov A. V., 

Wechsler R. H., 2005, ApJ, 624, 505 
Zentner A. R., Kravtsov A. V., Gnedin O. Y., Klypin A. A., 

2005, ApJ, 629, 219 



